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The  stress  function  for  a semi-infinite  elastic  strip  with  free  edges  is 
expanded  in  eigenfunctions  of  Papkovitch-Fad' le  type,  and  biorthogonal  func- 
tions are  constructed  by  use  of  the  adjoint  operator.  To  establish  complete- 
ness, the  same  expansion  is  obtained  by  a Fourier  transform  solution  of  the 
biharr  c n ic  equation . j)  £ - 0 0 1 f . 

When  end  data  is  supplied  in  the  form  of  prescribed  tractions  or  prescribed 


displacements,  the  coefficients  in  the  expansion  must  be  found  by  truncation  of 
an  infinite  set  of  linear  equations.  It  is  noted  that  the  methods  for  formulating 
such  equations  that  have  been  proposed  in  the  past  are  unstable  with  respect  to 
the  order  of  truncation. 

A systematic  search  among  all  possible  weighting  functions  in  the  biorthogonal 
family  for  use  in  a Galerkin  method  leads  to  a unique  choice  of  optimal  weighting 
functions  securing  maximum  stability  by  generating  a diagonally-dominated  matrix. 
Calculations  are  presented  for  a number  of  test  cases  using  this  method  together 
with  those  of  Benthem  (1963)  and  Johnson  and  Little  (1965),  and  with  a method  of 
direct  coll -cation.  The  method  of  optimal  weighting  functions  produces  consider- 
ably more  stable  values  for  coefficients  as  the  order  of  truncation  is  changed. 

This  opens  the  way  for  an  examination  of  the  convergence  of  the  eigenfunction 
expansions  for  data  of  this  class,  and  should  also  prove  useful  for  problems 
involving  cracks. 
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Key  Words:  Elastic  end  strip  problem,  Biorthogonal  eigenfunction 
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SIGNIFICANCE  AND  EXPLANATION 

The  problem  of  calculating  stresses  in  an  elastic  strip  |x^ | < 1,  x^  >0 
' from  given  end  conditions  has  proved  difficult  despite  its  practical  importance. 

The  Airy  stress  function  fi  possesses  an  expansion  in  eigen  functions,  but  the 
coefficients  in  this  expansion  can  be  found  explicitly  only  in  certain  particular 
cases.  When  the  end  data  takes  the  form  in  which  the  tractions  ^’12  are 

prescribed,  the  coefficients  have  to  be  found  from  the  truncated  form  of  an 
infinite  set  of  linear  equations.  These  equations  can  be  formulated  in  many  ways, 
but  it  has  been  found  that  those  proposed  in  the  past  lead  to  matrices  which  fail 
to  satisfy  the  criteria  for  stability  of  the  solution  as  the  order  of  truncation 
increases. 

A set  of  optimal  weighting  functions  for  use  in  a Galerkin  method  of  formu- 
lating the  linear  equations  has  been  found  by  systematically  searching  among 
possible  members  of  the  family  of  eigenfunctions.  These  functions  secure  maxi- 
mum stability  by  generating  a diagonally-dominated  matrix.  Calculations  are  pre- 
sented for  a number  of  test  cases  using  this  method  together  with  those  of 
Benthem  (1963)  and  Johnson  and  Little  (1965),  and  with  a method  of  direct  collo- 
cation. The  method  of  optimal  weighting  functions  produces  considerably  more 
stable  values  for  coefficients  as  the  order  of  truncation  is  changed.  This  opens 
the  way  for  an  examination  of  the  convergence  of  the  eigenfunction  expansions 
for  data  of  this  class,  and  should  also  prove  useful  for  problems  involving 
cracks. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC , and  not  with  the  author  of  this  report. 


MIXED  BOUNDARY  VALUE  PROBLEMS 
FOR  THE  ELASTIC  STRIP:  THE  EIGENFUNCTION  EXPANSION 


D.  A.  Spence 

1 . Introduction 

The  stress  function  ft(x,y)  in  a semi- in finite  elastic  strip  y > 0 

2 

with  free  edges  x = ±1  satisfies  the  biharmonic  equation  42=0,  and  pos- 
sesses an  eigenfunction  expansion  of  the  form 

_ -V 

n = l ck  e K wk(x)  (1) 

from  which  the  stress  components  are  found  as  the  second  derivatives  a , 

n , Si  . The  X.  are  the  zeros  with  positive  real  parts  of 
xy  yy  x 

X ± sinX  cosX  [ eigenfunctions j . 

Such  expansions,  associated  with  the  names  of  Papkovitch  (1940)  and  Fad'le 
(1940),  are  discussed  by  Lur'e  (1964),  Buchwald  (1964),  Buchwald  & Doran  (1965) 
and  many  other  authors. 

For  certain  combinations  of  boundary  data  on  the  edge  y = 0,  namely  when 
either  (i)  and  fl  , or  (ii)  and  are  both  given  on  this 

boundary,  the  coefficients  c^  can  be  written  down  as  quadratures  of  the 
boundary  data,  using  certain  functions  biorthogonal  to  appropriate  derivatives 
of  the  . The  first  of  these  combinations  was  treated  by  Smith  (1952); 

recently  Joseph  (1977)  has  discussed  the  convergence  of  the  resulting  expansion 
in  terms  of  the  boundary  data  for  this  case. 

In  the  present  paper,  we  consider  combinations  of  edge  data  for  which 
explicit  solutions  are  not  available.  First,  however,  expressions  for  the 
displacement  and  stresses  in  terms  of  Si  are  given  in  section  2,  and  the 
eigenfunction  expansion  (1.1)  is  developed  in  section  3,  the  stated  biorthog- 
onality properties  of  the  functions  u>k  .and  their  derivatives  being  derived 
formally  from  properties  of  the  differential  operator  associated  with  the  bi- 
harmonic equation  and  of  its  adjoint  operator.  This  derivation  does  not  prove 
that  the  expansion  (1.1)  is  complete.  However  the  solution  of  the  biharmonic 
equation  for  the  strip  is  again  obtained  in  section  4 by  use  of  the  Fourier 
sine  or  cosine  transform,  as  appropriate  for  the  boundary  conditions.  It  is 


found  that  the  transform  of  ft  (even)  has  poles  at  the  zeros  of  A + sin\  cosA, 
and  at  no  other  point.  The  residues  at  these  points  produce  the  expansion 
(1.1),  and  its  completeness  for  y > 0 can  thereby  be  established. 


In  the  practically-important  case  when  both  tractions  (i.e.  and  -n  ) 


are  prescribed  on  the  edge  y = 0,  an  explicit  solution  is  not  possible.  From 
(1.1)  we  have  the  expansions 


( nxyV=  0 

£ Ck 

(1.2a) 

^Uxx'y=0 

I Ck 

.(2)  • 

I’k  (x)  ; 

(1.2b) 

k ‘x21 

(x)  = 

( x),  but  the 

c,  must  now  be  found 
k 

where  <t>^  (x)  = 

by  solving  a truncated  set  of  equations  derived  from  these.  Various  ways  of 
deriving  such  a set  have  been  suggested:  Benthem  (1964)  equates  Fourier  sine 
and  cosine  coefficients  of  the  two  sides  of  (1.2a)  and  (1.2b)  respectively; 
Gaydon  and  Shepherd  (1964)  also  form  Fourier  coefficients  of  the  two  sides, 


but  evaluate  the  c^  by  optimizing  the  fit  to  a larger  number  of  coeffi- 


cients in  least  squares  approximation. 

A different  approach  involving  eigenfunctions  of  the  adjoint  operator  is 


due  to  Johnson  and  Little  (1965) . In  effect  they  use  certain  of  these  eigen- 


functions, say  i|i. 


(1) 


'Ia 


(2) 


Tm  ' Tm 

obtaining  an  infinite  set  of  equations  of  the  form 


, as  weighting  functions  in  a Galerkin  method,  thus 


m 


l . c,  + d 
mk  k m 


where 


c-  = [ F. 

jlUU>  .,(1)  x .(2)  .,.(2) 


(1.3) 


mk 


i|i"')dx 

vm 


(m^k)  , 


F =0,  and  d is  known  from  the  boundary  data, 
mm  m 


Johnson  and  Little  are  led  naturally  to  their  functions  ^ , i|>J^ 


by 


considering  the  biorthogonality  properties  of  a larger  set  of  functions 

1,2, 3, 4 which  gives  rise  to  the  explicit  solutions  mentioned 


(i),,li),  i 


k ' Tm 

earlier,  but  are  not  able  to  obtain  an  explicit  solution  for  the  present  case. 


and  it  can  be  shown  that  there  is  no  possible  choice  of  weighting  functions 


for  which  F , = 0 for  every  m^k , as  would  be  necessary  for  an  explicit 


solution. 
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J 


* 

J 


(M  \ 

To  ensure  that  the  solution  £ of  a truncated  set  of  2N  of  the  above 
equations  should  approach  a definite  limit,  the  solution  vector  of  the  infinite 
system,  it  is  necessary  for  the  matrix  to  be  diagonally-dominated,  in  the 
sense  that 

l If  j < | 1-f  | (1.4) 

k*n‘  mk  lmn 

for  each  m . A detailed  examination  in  section  5 shows  that  this  require- 
ment is  not  satisfied  by  the  matrices  arising  in  either  the  Benthem  or  the 
Jov  ison-Little  methods. 


Accordingly  we  seek  weighting  functions  that  will  at  least  comply  with 
(1.4).  For  this  purpose  we  examine  the  matrix  F when  iji  ^ , <|/2^  are  re- 


in m 


placed  by  functions  of  the  form 


X(1)  = A6(1)  + B*(3) 


X(2>  = C*(2)  + D*(4) 

m ym  *m 


It  turns  out  that  the  criterion  (1.4)  can  be  satisfied  with  just  one  choice 


of  the  four  constants  A to  D . The  weighting  functions  so  found  are 
referred  to  as  "optimal". 

The  methods  mentioned,  and  in  addition  direct  collocation,  i.e.  satisfying 

the  equations  (1.2)  at  a particular  set  of  points  on  (0,1),  have  been  tested 

by  application  to  three  representative  sets  of  initial  data  for  ft  and 

xx 

ft  on  y = 0;  the  results  are  compared  and  discussed  in  section  6. 
xy 
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2 . Stress  function  and  displacements 

To  define  the  problem  we  use  the  standard  suffix  notation.  Elsewhere  in 
the  paper  x,y  will  replace  x1»x2  as  coordinates. 

Consider  as  in  figure  1,  a semi-infinite  elastic  beam  occupying  the  strip 
“1  < xj_  < 1,  x2  > 0 of  the  x^  ,x2  plane. 


The  edges  >^-±1  are  treated  as  stress-free  throughout,  and  we  shall 
consider  only  distributions  of  stress  such  that 

°ll'a 22  are  even 


a^2  is  odd 


with  respect  to  • 


}(2.1) 


The  traction  on  the  end  x^ = 0 is  assumed  to  be  equilibrated  by  a force  2n| 
in  the  x 2 direction  at  infinity,  where 

N = t (a22)x2=°  dXl‘ 

The  stresses  may  then  be  written  in  terms  of  a stress  function  o satisfying 
the  biharmonic  equation  in  the  form 


°11  = n'22 

We  shall  write 

and  the  biharmonic  equation  is 
where  A is  the  Laplacian 


3x1 2,  3X? 


°12  _n'12'°22 


o + n 
li 


Att  = P 

AP  = 0 


(2.2) 

(2.3) 


1 Jf-2 

The  displacement  gradients  in  plane  strain  are  given  by 

2yUl,l=^  ” v^°ll  ” vo22'  2uu2  2=_V011  + ^ ~ v)o22 

p(ul,2  + u2,l>  = °12 
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(2.4) 


To  treat  boundary  value  problems,  it  is  convenient  to  introduce  in  addition  to 
P,  its  harmonic  conjugate  Q linked  by  the  Cauchy-Riemann  equations 


P,l  = 0,2  P,2  = -Q,l  <2'5> 

Q may  be  defined  as  zero  on  x^  = 0,  and  is  then  an  odd  function  of  . The 
shear  strains  can  then  be  written 


= -a  1?  * (1  - v)  q , 


The  vaJues  of  c^2,  o22  - N,  Q and  p on  the  boundary  Xj=  0 will  be  de- 
noted by  a 4-vector  f(xj)  with  components 


0 12 ' °22 


a22  - N 


f1  (AQ), 2 dXj^ 


■f(1)(x^ 
f (2)  (x,) 

f (3)  (X,) 
-f  (4)  (*1 1 


= f(Xj_)  (2.7) 


f^*^,  f^  are  odd  functions  of  Xj^ , f*2^,  f^4^  are  even. 


In  terms  of  these  functions,  the  gradients  of  displacement  along  the  edge  are 

/u  \ ( -f(2)  + (1-v)  f (4)  - 

- 1 


The  corresponding  expressions  for  the  case  of  plane  stress  are  obtained  by 
writing  yyy  in  place  of  v in  the  above.  These  are  equivalent  to  the  expres- 
sions (2.11  d)  of  Johnson  and  Little. 

The  biharmonic  equation  can  be  solved  explicitly  when  either  (a)  f ^ 


and  f'  ’ are  given  (i.e. 


and  u2)  or  (b)  f 


and  f 


are  given 


(i.e.  ct^2  and  u^) . We  proceed  to  find  the  solutions  for  case  (a) and  case  (b) 
in  turn,  by  use  of  Fourier  cosine  and  sine  transforms  respectively,  and  then 
combine  the  solutions  to  treat  (c)  (the  stress  problem) : o^2  and  °22  given, 

(d)  (displacement  problem)  u^  and  u2  given,  for  both  of  which  a unique  solu- 
tion also  exists,  but  can  be  obtained  only  in  terms  of  an  infinite  set  of 
linear  equations. 


3 . The  eigen  function  expansion 

From  now  on  we  shall  replace  x1(x2  by  x,y,  and  suffixes  will  be  used  to 
denote  derivatives. 

2 

The  biharmonic  equation  A S2  = 0 in  y>0  with  stress- free  boundary 


conditions 


= 0,  a = 0 on  the  edges  x = ±1 


and  symmetry  about  x = 0 possesses  eigen-solutions  of  the  form 
e w(x) 

(Re  A > 0) , where  w satisfies  the  fourth  order  ordinary  differential  equation 


2 2 2 

(D  + A ) u = 0 


> 


with  boundary  conditions 


x = 0: 


These  are  satisfied  by 


Du  = D (i)  = 0 


= Du  - 0 . 


in  = ( A cosA)  (tanA  cos  Ax  - x sin  Ax) 


when  A satisfies  the  equation 


A + sin  A cosA  = 0 


In  addition  to  the  root  A = 0 , the  equation  possesses  a conjugate  pair  of 


roots  A , A 
n n 


in  each  interval  (n  - i)  n < Re  A < (n-j-)ir 


V 


are  also 


roots.  The  asymptotic  location  of  the  roots  for  large  n is 

i in£  i i 

A , A • ^ i i i in  f + o(-4-  ) (3.5) 

n n 4 n £n  2 n 

where  £ = (4n-l)n.  The  first  10  roots  are  listed  in  Table  1. 

n 

Accordingly  we  seek  a solution  of  the  biharmonic  equation  subject  to  (3.1)  in 


the  form 


_ -*ky 

0(x,y)  = l c.  e 
k K 


“k(x) 


(3.6a) 


where  “k  (x)  = w(x, A^)  , and  the  c^  are  coefficients  to  be  determined  from  the 
remaining  boundary  data  on  x = 0.  In  this  summation,  A^  ranges  over  all  the 
roots  of  (3.4)  with  Re  A > 0 , i.e.  the  sum  means 

fThe  factor  (AcosA)-1  in  the  expression  for  w has  been  introduced  to  simplify 
later  expressions,  and  ensures  that  the  L,  norm  ||  u||  ■*  0 as  a ->■  °° . 
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Table  1.  Roots  of  (3.4)  in  first  quadrant 


n 

Re  A 

Im  A 

n 

n 

1 

2.106196 

1.125364 

2 

5.356269 

1.551574 

3 

8.536683 

1.775544 

4 

11.699178 

1.929405 

5 

14.854060 

2.046853 

6 

18.004933 

2.141891 

7 

21.153413 

2.221723 

8 

24.300342 

2.290552 

9 

27.446203 

2.351048 

10 

30.591295 

2.405013 

1 <lnen 

The  asymptotic  expression  j £n  - — - gives  6 figure  agreement  with 

n 

1 2 2 

Re  A at  n = 4,  and  ■=-  in?  + [{£nt  ) - 2s.n5  -l]/5  agrees  with  Im  A to 

n £■  n n n n n 

4 x 10-5  at  n = 9. 


I 
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(3.6b) 


We  find  that  L = L1  , with  boundary  conditions 
y = o:  ^ = *<3>  = 0 

y = l:  *(3)  = 0,  D(*(1>  - *<3>> 

These  are  satisfied  by 


II 

•e- 

OJ 

II 

o 

=0,  D(^(1) 

- <l> 

(1)  \ / 

A (3) 

I - 1 [ 

(3)  )'  2 V 

4(1) 

(3.11) 


\ *(3)  / 2 V *(1)  J , (3'12) 

where  A is  a root  of  (3.4).  The  normalizing  factor  j has  been  introduced 


so  that  the  scalar  product 


U<1,3)  (1,3)  _ a . 

m n mn 


(3.13) 


(2  4 ) 

The  vector  $ ' is  also  an  eigenfunction  of  the  same  matrix  operator 


with  the  different  boundary  conditions 


y = 0:  D<|>(2)  = D<f> ( 4 3 


y = 1: 


(2). (4) 


DU(2)-  <t>(4)  ) = 0. 


(3.14) 


In  this  case  the  adjoint  is  again  L , with  side  conditions 
y = 0:  Dii(2)  = D i|<(4)  =0 


y = 1: 


<i>(4)  =0,  D \|>(43  = 0 


whence 


and  likewise 


*(4) 

(2)  , (4 ) 


(2,4)  *(2,4).  = 


) = « 


(3.15) 


(3.16) 


The  biorthogonality  relations  can  be  confirmed  by  direct  calculation  (Appendix  a)  . 


Writing  < , i = 1,2 ,3, 4 we  find  that 

m n mn 


(3.17a) 


while  for  m^n 


„<1>  - -v<3>  - 2X"Xn  v(2)  _ „(4) 

mn  mn  ^z  ^z  mn  mn 

m n 


2 A2  c 

m mn 

2 2 

A - A 

m n 


(3.17b) 


c = (cos  A cos  A„)  ( cos  A x cos  Ax) 

ran  m n nv  n 


(A  tan  A - A tan  A )A\2  - A2) 
m m n n/  m n 


(3.17c) 


2 -1 

From  the  equation  tan  A + A tan*  +1=0,  which  is  a form  of  (3.4) , we 
deduce  that 

tan  A = ii  (1  - 1/4A2)1/2  - 1/2 A , 
m m m 

the  upper  and  lower  signs  corresponding  to  Im  A < 0.  Thus  tan  A^  -*•  ±i  in 
upper  and  lower  half  planes,  and 

^ ciwJ  % n(m+n)  ' ir(m-n) 

accordingly  as  the  imaginary  parts  of  *m'*n  have  the  same  or  opposite  signs. 
Explicit  solutions 

If  the  functions  f^^  and  are  prescribed,  the  biorthogonality  relation  (3.13) 


immediately  gives  the  coefficients  in  the  expansion  (3.6)  as 


ck  ’ < 4"  ^ 

+ l/-j[3)  f (3)  > . 

(3.18a) 

Likewise  when  f^)  and  f^4^ 

are  prescribed 

„ , , (2)  ,(2) 
Ck  = ( \ f 

+ ^4>  f<4>>  . 

(3.18b) 

In  other  cases,  biorthogonality  does  not  lead  to  closed  expressions  for  the 
(Cj^ } , although  the  (i|;^^  } provide  useful  candidates  for  weighting  functions 
for  use  in  a Galerkin  method.  We  return  to  these  problems  in  Section  5,  but 
before  this  in  section  4 we  demonstrate  the  completeness  of  the  expansion  (3.6) 
by  use  of  Fourier  transformation  methods. 

Eigenvectors  of  the  Fredholm  determinants  ^ 

We  conclude  this  section  by  noting  a set  of  identities  among  the  coeffi- 
cients K ^ <t>^>  that  can  be  inferred  as  follows  from  the  expan- 

mn  vm  *n  r 

sions  (3.7)  and  the  explicit  solutions  just  obtained.  In  the  case  i = 1 , we 
have 


f(1,(x)  - f1(  l *11)  (x)^1)t)f(1)  (t)dt  = /15V1)  (x)^3)(t)f  (3,(t)dt.  (3.19) 

« m in  a iu  in 

U . U (1) 

This  is  in  the  form  of  a second  kind  Fredholm  equation  for  f (x)  , given 

f(3).  However,  f(1)  and  *(3) 


f represent  independent  boundary  information. 

It  is  physically  obvious  that  neither  can  imply  the  other.  By  the  Fredholm 
alternative  theorem  however,  the  equation  (3.9)  has  a solution  for  every 
f^3'  unless  the  homogeneous  equation  obtained  by  setting  f^3^  =0  possesses 
a non-trivial  solution.  We  therefore  expect  such  an  eigen-solution  to  exist. 


-10- 


(3.20) 


.21) 


(3.22) 


.23) 
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Fourier  transform  solution 


In  this  section  we  seek  the  solution  of  the  biharmonic  equation  in  the 

strip  |x | _<  1,  y > 0 by  use  of  a Fourier  transform,  subject  to  the  condition 

n = 0 , n = 0 on  x = i 1 (4.1) 

xy 

and  given  data  on  the  end  y = 0 in  either  of  the  forms 

(a)  SI  = f (1)  (x)  , /X(Afl)  vdx  = f (3)  (x)  (4.2) 

xy  0 y 

or 

(b)  fixx  = f (2) (X) ' *n=f(4)(x)  (4.3) 

Only  solutions  for  Si  even  in  x will  be  treated,  and  we  assume  that  the 

boundary  data  is  consistent  with  this,  i.e.  that  f*3^  and  f^  are  odd 

with  respect  to  x , and  f*2^  and  f ^ even,  and  further  that 

/ f'  *(x)dx  = 0,  in  accordance  with  (2.6). 

0 


We  carry  out  the  solution  in  turn  in  the  two  cases. 

Case  (a) : Multiply  the  biharmonic  equation  on  both  sides  by  cos  £y  and 
integrate  with  respect  to  y from  0 to  °°,  obtaining 

(D?  - 52)2  n(c)  = <flyyy  + 2flxxy  - c\)y=o  (4 

where 


Si  (c)  (x;£  ) 


= J°°  (cos^y)  n (x,y)d 

0 


* 


and  D = 


,(c) 


W 


as  before.  The  solution  of  this  equation  even  in 


---hi 

2C  0 


X k(1'3) (Cx-et) f (1,3) (t)dt  + f1  f 


(1) 


(t)dt 


(4.5) 


+ A ( O coshCx  + B ( Ox  sinh^x 
where  k^3*  (s)  = -2cosh  s,  k^3^  (s)  = -s  sinh  s, 

and  k ' 3>  f (1 ' 3)  denotes  the  scalar  product  k^'f^  + k^f*3^  ' The 

functions  A,  B are  to  be  determined  from  the  conditions  (4.1),  which  can 
be  applied  in  the  form 

a (c)  = 0 , fi(c)=  0 on  x = 1 

X 
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giving 


M(?)  ( *)  = -V  / 


1 ( k(1,3)  (f-r.t) 


2?2  0 


where 


k'  (1,3>  (?-? t) 


If (1,3) (t) dt 


(4.6) 


M(  ?) 


cosh  6 sinh  ? 


-1 


sinh  ? cosh?+?  sinh  6 


Thus 


(c) 


is  given  by  (4.5)  , with  the  last  two  terms  expressible  as 


A cosh  ?x  + B*  sinh  ? x = 


k(1,3) (?-?t) 


(4.7) 


Jl( cosh?x,  x sinhCx)  M-1  (?)(  .. 

2?z  0 V*’  (?-?t) 

The  inverse  matrix  in  the  integrand  is 


f <1,3) (t) dt , 


, . f cosh?  + ? 3 sinh?  -sinh? 

M (?)  = (det  M)  1 


-sinh? 


cosh? 


detM=  1 + (1/?)  cosh?  sinh?. 

a'~'  is  an  even  function  of  ?,  and  a detailed  examination  shows  that  it  is 
regular  at  ? = 0.  The  inverse  cosine  transform  is  therefore 


with 

, (c) 


0(x,y)  = — l°°  a (c)  (costy)d?  = i-  J”  nlc)e1^Yd? 


= I r n(c)  i?y. 


(4.8) 


(c)  -2 

As  ft  = 0(£  ).  For  future  use  the  asymptotic  form  is  quoted  in  full: 


fi(c)  v -±-  !l  [e_c ' x_t i (-2f (1)  (t)  + ?|x-t| f(3)  (t))sgn(t-x) 


4 ?2  0 


+ e"5(X+t) (-2f(1) (t)  + ?(x+t) f(3) (t) ) ] dt  (4.9) 


+ f1  f (1) (t) dt  . 

? x 


The  integral  (4 .8)  can  therefore  be  evaluated  for  y > 0 as  2i  times  the 


sum  of  residues  of  the  integrand  at  its  poles  in  Im  ? > 0 . The  only  such 

-1 


poles  are  those  of  M (?)  at  the  zeros  of  det  M(?),  i.e.  at  ? = iA  , where 
COO  = 2>  + sin  2 X = 0 , Re  A > 0 (4.10) 


Each  of  the  eigenvalues  A , given  in  (3.5)  produces  a contribution 
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- A 


n y j ^ 

( / (cos  Ax , i x sinAx)  ( resf 

i.\  0 


M_1>  ( i.‘a:-nt))^t)dt 


The  residue  of  M can  be  expressed  as  the  matrix  product 


(res  M 3) 


iA 


( tan  A 1 1 , i. 

I i J (-  2 tanA,-  j) 


5=iA  cosA 

and  the  integrand  (4.11)  then  factorizes  since 
, -1 , , . . , . rtan  A 


(4.11) 


(4.12) 


(AcosA)  (cosAx,-  ix  sinAx)  ( ^ ) = w(x,A) 


' k(1,3)  (iA-iAt) 
\k ' (1,3) (iX-iXt) 


where 


(-  i tan  A,-  ^-) 

00  = (AcosA)  1 (tanA  cosAx  - x sinAx) 

iJ/^(t,A)  = (cosA)  1 sin  At 

i^3^  ( t . A ) = j (cosA  sinAt  - t sinX  cosxt) 


= U(1)  ( t , X ) , * 


(3) 


(t , X)  ) 


(4.13) 


J 


are  precisely  the  functions  already  encountered  in  Section  3.  Thus  the  solu- 

-\ 

V 


tion  of  the  (1,3)  problem  in  y > 0 is 

Xk 


V "XkY 

ft  (x,y)  = I c ( A.  ) e (»)  (x , Av ) 


where 


(4.14) 


c(xk)  = (*(1'3) (t,xk) , f(1,3) (t)) 

A,  in  the  summation  ranges  over  A ,A  , n = 1,2,...  and  the  oo(y,X,  ) 
k n n K 

form  a complete  set  of  eigen  functions  on  (0,1). 

Moreover  the  expansion  (4.14 (holds  on  y = 0.  This  can  be  seen  since, 
for  y < 0,  the  inverse  Fourier  transform  (4.8)  is  evaluated  as  ( — 2 i ) times 
the  sum  of  the  residues  in  Im  X < 0 . These  residues  are  the  negative  of  those 
in  Im  X > 0 . Hence  the  expression  (4.14)  holds  with  | y|  in  place  of  y , 
and  is  continuous  at  y=  0. 

(1  3) 

The  biorthogonality  already  established  between  the  i)<  ' and  the  odd 

derivatives  ^ 1 ' 3 3 follows  immediately  from  (4.14)  since  it  gives  the  expan- 
sions 


(1) 

" (fJxy)y=0  ~ l ck 

,(D 

♦k 

(x)  , 

*(1)  = 

(3) 

= (/*  (*ft)  ydx)yto=  ! 

l ck 

4 (x), 

*(3)  = 

= (A_1D3+XD)u 


(4.15) 
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From  (4.14)  we  deduce  that  ( * j 1 ' 3 * , ' 3 ' ) - • 

In  case  (b)  the  steps  are  similar:  given  f 

sine  transform  of  the  biharmonic  equation,  obtaining 


and  f 


(4) 


(D2-e2)2o‘s)  = u3n  - 2enxx  - en  )y.0 


n 


where 

with  solution 
ft(s)  = - 


(4.16) 

we  form  the 


(4.17) 


(S) 


(x;£)  = /“ (sin^y) ft (x,y) dx 
0 


1 - k(2,4)  (Ex-Ct) f (2,4)  (t)dt  + 7 /1max(x,t) f (2)  (t)dt 

2C2  0 5 0 

+ Acosh  5x+  Bx  sinh  £x 


(4.18) 


where  k(2)(s)  = 2 sinh  s,  k(4)(s)  = s cosh  s - sinh  s.  The  term 
Acosh^x  + Bx  sinhC  x is  given  in  terms  of  f^2,4^  by  an  expression  identical  in 


form  to  (4.7).  In  this  case  ft 


(s) 


reference  we  quote  the  asymptotic  form  of  ft 


is  odd  with  respect  to  £ 
(s) 


For  future 


as 


!(s)  v ~ J1[e  (2f(2)-(l+s|x-t|)f(4)) 

+ e"5(x+t)  (2f(2)  - (l+C(x+t))f(4))]dt 


(4.19) 


When  steps  identical  to  those  for  ft 
(4.14)  is  again  obtained  but  now 


(c) 


+ j /1(t-x)f(2)  (t)dt 

X 


are  applied  to  (4.18)  the  expansion 


,,(2,4)  f (2 , 4) 
ck  = (*k  ' f ’ 


,(2)  y _ .U)  - V _ 

f = 1 ck  ' f ' iCkTk 

the  biorthogonal  functions  in  this  case  being 


(4) 


(4) 


<l> 


^2*(t,X)  = - (1/cosx) cosxt. 


*(4)  (t,  X)  = -i-X2  01  (t,  X) 


and  the  second  biorthogonality  relationship 


(*; 


is  again  obtained. 


<2'4),  *<2'4))  = « 


ik 


(4.20) 

(4.21) 

(4.22) 


5 . Mixed  boundary  value  problems 

In  this  section  we  consider  the  problem  of  finding  the  coefficients  {cn } 

in  the  expansion  of  n and  its  derived  functions  when  the  data  on  the  boundary 

does  not  fall  into  one  of  the  types  ((f^,f^)  or  (f^,f^4')  given)  for 

which  an  explicit  solution  was  found  in  the  last  section. 

Two  cases  are  of  particular  practical  importance: 

(c)  The  stress  problem:  and  o^2  are  prescribed  on  y = 0,  i.e. 

f^(x)(odd),  f ^ (x)  (even)  given.  (5.1) 

(It  will  be  assumed  that  f(1)(0)  = f(1)(l)  = 0,  and  f1  f(2)(x)dx=  0 , to 

0 


conform  with  (2.1).) 

(d)  The  displacement  problem  in  which  u and  v are  prescribed  on  y = 0, 
i.e.  the  combinations 

t *3  \ / O \ / O \ f A \ 

(5.2) 


gd)  _ f(l)  + (1_v)f(3)  (odd),  g(2J 


f(2)  - (l-v)f(4>  (even) 


are  given  . 

We  note  first  that  explicit  solutions  to  these  problems  cannot  be  found, 
i.e.  in  case  (a)  there  is  no  set  of  functions  biorthogonal  to  the  set 


(1)  .<«> 


(-A  Du>  / D (i)  ) . 
n n n 


If  there  were  such  functions  (Z ^ ,Z ^ ) , we  should  have  for  each  m ^n . 

m ' m 


(Z<1'2),*<1'2))  = 
m n 


: ( X Z (1) + DZ (2) ) Du  > 


(on  integration  by  parts) . This  can  hold  for  every  n only  in  the  trivial 


.(1)  . 
Jm 


(2) 

Z = constant, 
m 


The  coefficients  for  case  (c)  can  therefore  be  found  in  principle  only 
by  deriving  an  infinite  set  of  linear  equations  from  the  expansions 


f (1)  (x) 


: ( 2 ) 


I 


(1) 


(X) 


(X)  = 5 cn  *n 


(2) 


(x) 


(5.3a) 

(5.3b) 


and  hoping  to  recover  the  cn's  with  sufficient  accuracy  from  a truncated 
set.  In  these  equations,  the  summation  ranges  over  each  eigenvalue  Xn  in  the 
first  quadrant,  and  its  conjugate  X . The  corresponding  coefficients  and 
eigenfunctions  are  also  conjugate. 

Three  essentially  different  methods  for  formulating  equations  have  been 


examined . 


-16- 


r 


(A)  Direct  collocation:  In  which  a truncated  set  of  the  equations  (5.3  a,b) 
are  identically  satisfied  at  a suitable  set  of  points  (x  } of  the  interval 
(0,1)  . 

(B)  Galerkin  methods  Here,  sets  of  weighting  functions  W*3^,w'2'  are 
m m 

introduced  and  mean  equations 

<’^(1))  = I^’^n^n'  ^2)'  f'2))  = K^^n  <5’4> 

are  obtained  and  solved  simultaneously.  The  simplest  such  equations  are  pro- 
vided by  the  set 

t ■>  \ 

(5.5a) 

These  were  proposed  by  Benthem  (1963),  and  are  equivalent  to  those  used  by 


= sin  mux,  = cosimrx. 

m m 


Gaydon  & Shepherd.  With  our  notation  the  matrix  elements  are  given  by 


< sin  m^x 


(1) 


) = < 


cos  mirx 


2 

. 2mnX  tanX 

( 2 )._,..  m- 1 n n 


> = (-1) 


.71 — TT71 

(X  -mu  ) 
n 


(5.5b) 


(C)  Use  of  iorthogonal  functions 

Johnson  & Little  in  effect  add  equations (4 . 14 ) and  (4.20)  to  obtain 


2c 


m 


(^1'2),f(1’2))  + (^3'4),f<3'4>). 


(5.6) 


Here  f^3'2^  is  known,  but  f^3'4^  must  be  calculated  as  \ c^t^3'4) 
From  (3.17) 

f 1 (m  = n)  r 1 


(4l3'4),*i3’4))= 


rr 


n 


-2X  (X  tanX  - X tan  X ) 
m m m n n__ 

(X  - X ) 2 (X  +X  ) 
m n m n 


(m  ? n)  3 


(5.7) 


say.  This  gives  the  set  of  equations 


where 


c = y F c + d 
m *■  mn  n m 


d = (*<3'2>,f<3'2> 


(5.8) 


m rm 

It  is  to  be  noted  however  that  this  set  of  equations  can  also  be  obtained 
directly  by  a Galerkin  method  from  (5.3a)  and  (5.3b),  using  the  weights 
W(1'2)  = u,(1'2)  and  adding  . 
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For  an  infinite  set  of  equations  of  the  form  (5.8),  Kantorovich  & Krylov  shew 


that  the  solution  (cn)  is  stable  in  the  sense  of  approaching  a limit  as  the 

number  N of  equations  become  large  if  the  norm  ||  F|  1^5  sup  £ |f|  is  - 1 . 

n 

However  this  condition  is  not  satisfied  by  the  coefficients  (5.7);  in  fact  the 

sum  I If  | tends  to  infinity  with  m . 

‘ 1 mn1  2 

n 

To  discuss  the  sum,  it  is  convenient  to  write  the  equations  more  precisely 


7 (F  c + F*  c ) + d 
L , mn  n mn  n it 

n=l 


where 


F*  = -2  A (A  tanA  -A  tanl  )/(A  -A  ) 2 ( A +A  ) . 

mn  mm  mn  n'  mn  mn  J 

The  summation  now  extends  over  X^  in  the  first  quadrant  only,  and  we  need 
only  consider  X^  in  the  first  quadrant  since  the  conjugate  equation  holds  for 
Xm  . The  asymptotic  expressions  (3.5)  show  that 


X ,X  v mit,nu,  while  tanX  , tan  X -*•  i. 
m n m n 


Therefore 


F v - 
1 mn 1 it 


I t*  I 2m  , . —2 

lFmnl  ' T (m-n) 


(m  ? n) 


(5.10) 


(For  m = n,  F =0,  F ^ 2imw/  (in4mtt)  ) . The  imaginary  part  of  F does 
mn  mm  mm 

not  appear  in  the  diagonal  when  the  equations  (5.8)  are  formulated  in  real 


and  imaginary  form.  Then 


y If  | 't/  — in  2m,  y |F  I 

L ' mn 1 it  ' L 1 mn 1 


m fu  1 . x - mn 

T 11+  + ~ — * 


Similar  behaviour  is  exhibited  by  the  elements  (5.5b)  of  the  matrix  formed 


by  sine  and  cosine  weighting  functions,  and  we  might  expect  this  to  be  re- 
flected in  instability  of  the  solution  of  truncated  sets  of  these  equations. 

In  fact  this  was  shown  in  numerical  experiments  described  in  the  next  section. 

2 

It  is  evident  that  the  difficulties  arise  from  the  factor  (X  -X  ) m 

m n 

the  denominator  of  F . 


D.  Optimal  weighting  functions 

In  the  hope  of  producing  a diagonally-dominated  matrix,  we  may  consider 
weighting  functions  of  the  general  form 


(5.11) 


This  form  is  a natural  extension  of  the  biorthogonal  weighting  functions 

j i^2'etc.  anc*  contains  sufficient  disposable  constants  to  ensure  the 

2 

suppression  of  the  factor  (X  -X  ) in  the  denominator  of  the  matrix  elements. 

m n 


(1) 

m 

= A $(1)  + 
m 

B 

♦ <3) 
m 

(2) 

m 

= C *(2>+ 
m 

D 

*(4) 

m 

If  we  write 

f(1)  + af<3)  . g(l)f  f(2)  _af(4)  = g(2)  (5.12) 

then  the  stress  (u=0)  and  displacement  (a=l-v)  problems  may  be  treated  to- 
gether as  cases  in  which  the  functions  g^  , g1'2'1  are  known  on  0 < x < 1 . 
We  shall  therefore  seek  optimal  weighting  functions  in  the  above  sense  for 
g*1'  and  g*2*,  and  expect  the  result  to  depend  on  the  parameter  a.  These 


produce  the  equations 


T G c + d 
L mn  n m 


where  d = (x^'2^»  g^'2^)  is  known,  and 
m m 

- _ , , (1),  '1)  . . (3) . . (2 ) . ( 2)  . (4),. 

G = ~ < X_  (<l>  + 0 4 ) + x_  V — ot  4>_  )> 

mn  mn  in  n n m n n 

Results  in  appendix  A show  that  for  m ¥ n 


(5.13) 


4 X X 

m n 

(X2-X2 
m n 


A(X2  + X2 ) + 2CX  X 
m n m n 


- B + A + 


o 2 

a D + tXC 
m n 


where  c 

rr 

as  before  , and 


- 4 a (B-D)c  + 4aB  t 

mn  mn 

c = (X  tan  X - X tan  X )/(x2  - ) 
mn  m m n n m 


t = (tan  X + tanX  )/(X  + x ). 

mn  m n m n 

To  suppress  the  factor  (X  - X )3  in  the  denominator  we  must  have  A = -C. 

m n 

Without  loss  of  generality  we  can  choose  C = 1,  A = -1  when  the  expression 
reduces  to 


- 4 


A A 

m n 


<x  +a  r 
m n 


-^B+a ) A A + A2D  + A2  a 
. mm  m n 

+ ? 5 + ot(B-D) 


m 


n 


c ^ + 4 uB 
mn 


tanA  +tanA 
m r 


A + A 

m n 


and  if  B = | - a , D = - a 


the  remaining  factors  of  (a  -A  ) in  the 

m n 


denominator  are  suppressed  leaving 

Gmn  = 2VXm  tan  Xm  ' Xn  tan  Xn)/(Xm  + V 


+ 2a (l-2a) (tan  Am  + tan  A )/(A  + A ). 

m n m n 


(5.14) 


This  formula  holds  when 
the  diagonal  term  is  zero  when 


xm  * xn 


Further  results  in  the  appendix  show  that 
0,  and  otherwise  is 
( tanA 


Gmm  = 4 “ + 2a (l-2a) 


m / 


(5.15) 


The  weighting  functions  are  then 

-l1’  - 


(5.16) 


(2)  = J2) 


(|  + a) 


><4> 

Km 


For  the  case  a - 0,  i.e.  the  stress  problem,  case  (c) , the  diagonal  elements 
are  all  unity , and  the  sum  of  the  absolute  off  diagonal  elements  is  strictly 


less  than  1.  This  can  be  seen  as  follows:  Again,  denote  by 

of  (5.14)  with  Ar  in  place  of  A_.  Then  for 
* 


of 


mn 


mn 


n - - n* 

are  respectively 

2im  m-n 


a 


the  value 


0,  the  asymptotic  values 


whence 


(m+n) 


(m  + n)'2 


(5.17) 


m 


I <K 


I)  - 2-t  l 


2m 


n=l (m+n) 


+ — i 

IT  u 


2n 


n=m+l  (n,+n) 


= y 


3 i f- 


n=l  (m+n) 


l \ r 

n=m+l  1^  (m+n) 


* 1 

(m+n)  J 


as  m -*• 


3 2 

The  estimate  — is  improved  to  — when  the  real  and  imaginary  parts  of 
the  equation  are  separated. 


When  a is  different  from  zero,  it  is  possible  to  find  a positive  constant 
y depending  on  a such  that  £|  ( jy)YGmnl  < I 1— 4c» | ; this  implies  that  provided 


I m ^d 


-*  0 as  m -*•  °°,  the  equations  (5.13)  possess  a unique  bounded  inverse. 
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6.  Calculations 


6.1.  End  Stresses  Prescribed 

The  four  methods  described  in  the  previous  section: 

A Collocation 

B Sine  and  cosine  weighting  functions 
C Biorthogonality  relations 
D Optimal  weighting  functions 

were  programmed  and  applied  to  the  stress  problem  (c)  for  three  test  cases: 

- 3x2 


(1) 

f(D 

= 0 

, f(2)  = 

(2) 

f(1> 

= 0 

, f(2)  = 

(3) 

fU) 

= X 

- x3  , f 

(f(1)  = <°12 

^y  =0 ' 

f (2) 

= (o22  > y 

eigenvalues 

V- 

* XN' 

Xl* • ’ XN ’ 

were  separa 

ed  to 

produce  real  i 

H 

(2N*2N) 

where  now 

z2n- 

1 + 1 

z_  = c 
2n  n 

( 1 (0  <_  x < |) 

V-l  (i  < x < 1) 


(6.1) 


(2)  _ 


The  real  and  imaginary  parts  of  the  equation 


(6.2) 


Similarly  the  right  hand  sides  were  separated  as  dn 
method  C,  (5.9)  takes  the  form  (6.2)  with 


b_  .+  i b_  . 
2m- 1 2m 


M 


M 


2m- 1 2n-l 


Re  ( F + F ) 


2m- 1 2n 

f12m 
M 


= Im ( F-F  ) 


mn 


(6.3) 
Thus  in 


) (6.4) 


_ , = -Im(F+F  ) 

2n-l  mn 


Re (F-F  ) 


2m  2n  mn  mn  j 

N was  given  the  values  5,  10,  20  and  40  in  turn.  The  resulting  coefficients 

c^  are  listed  in  Table  2 for  cases  1 and  2.  Case  3 produces  results  qualita- 
tively very  similar  to  those  in  case  1.  The  coefficients  have  been  used  to 
calculate  the  distributions  of  normal  and  shear  stress  at  distances  c from 
the  boundaries,  using  the  augmented  expansions 
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mm 


r 


- \ y n > 

a12(x,y)  = l cn  e *n  (X) 

n 

° 22  (X,y)  = l °n  6 ^ *n2)(x) 

truncated  to  2N  terms.  The  results  of  these 
with  the  various  methods  are  indistinguishable 


calculations  with 
except  at  x = 0 , 


(6.5) 

N = 40 
and  are 


plotted  in  figures  3,  4,  5. 

At  y = 0,  the  accuracy  with  which  the  expansion 
data  was  assessed  by  computing  o,  where 


2 

o 


wi th  x . = 

l 


x~2 

M 


eluded  in  Table 


1 

M 


M 


l 

i=l 


<v 


(1) 

n 


for  M = 100.  The  resulting  values 

2. 


(5.3)  fits  the  input 

(2)<V  - l Vn2)(xi>)5 

of  o are  also  in- 


6 . 2 End  displacements  prescribed 

In  addition,  the  optimal  weightincr  function  approach  was  applied  to  the 
displacement  problem  for  a beam  bonded  to  a wall  with  u=0,v=0  on  x=0, 
0 < y < 1,  under  unit  tension  at  infinity.  For  these  boundary  conditions  (2.7) 
and  (2.8)  give 


g 

g 


(1) 

(2) 


+ 


(l-v)f 

(l-v)f 


(3) 

(4) 


0 


- v . 


(6.6) 


Wfe  denote  this  case  by  E4 , and  have  obtained  the  coefficients  when  v = 0.3. 
These  are  listed  in  Table  3. 


Discussion  of  results 

For  cases  1 and  3 in  which  the  data  is  smooth,  with  N = 40  methods  B and 
D produce  virtually  identical  results,  with  agreement  to  4 or  5 signifi- 
cant figures  for  the  earlier  coefficients  and  the  same  rate  of  decay.  Methods 
A also  agrees  to  3 or  4 figures  for  the  first  N/2  coefficients,  but  is  less 
accurate  beyond  this  point.  Method  C (Johnson  and  Little)  is  significantly 
less  accurate  than  either  of  these,  and  in  fact  becomes  worse  as  N is 

increased.  There  is  a significant  difference  in  the  rate  of  convergence  of 

(N) 

the  solution  vector  z with  methods  A,  B and  D:  for  method  D the 

first  10  coefficients  are  already  close  to  their  final  values  when  computed 
with  N = 10.  This  is  not  so  with  collocation,  method  A,  and  the  stability 
is  less  marked  with  the  fourier  coefficients  method  B.  The  biorthogonality 
equations  C are  again  significantly  worse  from  this  point  of  view. 
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The  corresponding  values  of  a are  of  the  same  order  of  magnitude  fop 

-2  -2  -3  -4 

given  N for  methods  A,  B and  D,  (i.e,  10  ,10  ,10  , 10  for  N = 

5,  10,  20,  40)  but  substantially  larger  in  the  case  of  method  C , and  in 

( 2 ) 

fact  increasing  with  N in  this  case.  For  case  2,  in  which  f is  dis- 
continuous at  y = 0.5,  the  coefficients  decrease  very  slowly.  The 

convergence  is  discussed  in  the  next  section.  The  optimal  weighting  func- 
tion method  D again  provides  the  most  stable  solution  as  n increases, 
and  the  Johnson-Little  method  C the  least  stable,  but  in  this  case  the  col- 
location method  A is  more  successful  than  method  B in  its  agreement  with 
the  earlier  coefficients,  for  which  we  may  assume  that  method  D gives  the 
best  estimate. 


With  all  four  methods  in  case  2 the  values  of  the  goodness-of-f it  para- 
meter o are  much  larger  than  in  case  1,  but  are  smaller  (0(10  ^))  with 
methods  A and  B than  with  method  D,  where  a is  0(1)  . For  method  C , 
o becomes  larger  as  N increases,  and  is  6.6  for  N = 40.  A print  out  of 
the  sums  of  the  series  (5.3)  using  the  calculated  coefficients  shows  that 
with  method  A the  fit  is  extremely  close  (3-4  figures  accuracy)  except 
in  the  immediate  neighbourhood  of  the  discontinuity,  whereas  with  method  D 
the  fit  is  poor.  This  reflects  the  fact  that  method  D obtains  the  coeffi- 
cients in  the  expansion  of  the  biharmonic  function  0 , but  that,  although 
this  expansion  is  convergent  at  y = 0 , the  expansions  (6.6)  of  the  stress 
functions  are  convergent  only  for  y > 0.  However  the  remaining  methods 
provide  interpolations  at  y = 0 , of  which  collocation  is  the  most  success- 
ful. 

Convergence 

A print  out  of  the  ratio  of  sums  of  absolute  off-diagonal  elements  to 
the  diagonal  element  in  each  row  was  made  for  method  D and  E,  and  confirms 
as  expected  from  theoretical  considerations  the  ratio 


l 

j»  j^i 


M . . I / I M • ■ I 

l]  1 1 'll  1 
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in  all  cases  decreases  steadily,  e.g.  from  about  .6  to  .2  for  the  stress 
problem  (2)  when  N = 40  (80  equations). 


The  right  hand  side  of  the  equations  for  method  D case  2 is  also  in- 
cluded in  table  2,  and  shows  that  c d as  n becomes  large,  because  of 

n n 

the  diagonal  dominance.  The  effect  would  be  much  more  marked  for  elements 

that  decay  more  rapidly,  and  it  appears  that  the  asymptotic  behaviour  of  the 

c can  be  estimated  from  that  of  the  d in  this  case.  Thus  the  biortho- 
n n 

goral  weighting  functions  provide  the  basis  for  a discussion  of  the  con- 
vergence of  series  (5.3).  In  the  three  cases,  the  expression  for  d^  for 
this  method  are 


x(1) 

f(1)^ 

Am 

/ 

tanX 

-2  1 

m 

X 

1 " X2 
m 


^ 0 (m  X) 


- 


l ) 

( O -L 

if  cos  2 

1 

4 

'v  m 

tanX  ^ 

\ | 

( ! _ 5 \ 

l X2 

\ in  1 

) 1 

\ Xm) 

(|  cos  i Xm  - i cos  | xm  - 1) 


18 
, 3 


(case  1) 


(case  2) 


(case  3) 


% 


0 (m  2) 


/tanX  \ 

(-v)  (l+2a)  ( % o (m  ) (case  E4) 

In  Appendix  B,  expressions  are  listed  for  the  L ^ norms  of  the  eigen- 
functions. The  asymptotic  behaviour  is 

II  ♦^1)||  - II  ^2>ll  - §n1,/(ien4n*)3/2  , 

while  ||  u>n  II  ^ (2nw)  1(jin4m)  3/'2  . 
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— w 


1 


— 2 5 

Hence  when  | | n,  n as  in  case  ( 1 ),  the  majorant  series 

nc„in,'1>ii 

is  convergent;  in  case  (2)  on  the  other  hand,  |c  I converaes  only  like  n °' 

n 

-a  y 

and  the  individual  terms  of  the  series  increase.  However  the  exponentials  e ri 
ensure  that  in  all  cases  the  series  (6.6)  is  convergent  for  x > 0 . 

Cor  elusions 

The  conclusions  to  be  drawn  from  these  results  are  that: 

(i)  The  optimal  weighting  functions  (5.16),  leading  to  the  equations 

(5.13)  with  matrix  G given  in  (5.14),  (5.15)  provide  the  most  stable 

system  of  equations  from  which  to  determine  the  coefficients  c 

n 

(ii)  The  use  of  unmodified  biorthogonal  weighting  functions  is  unsatis- 
factory, and  becomes  more  inaccurate  as  N is  increased. 

(iii)  In  cases  involving  discontinuities  in  edge  data,  collocation 
provides  a better  interpolation  than  Galerkin  methods. 
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Appendix  A : Quadratures  of  eigenfunctions 
The  quadratures  quoted  in  the  paper  are  derived  from  the  following  re- 
sults, all  of  which  are  obtained  by  integration  by  parts  and  use  of  the 
boundary  conditions  (3.2)  on  u and  its  derivatives,  and  of  the  fact  that 
X.  = -sin  X nos  1 ^ 3 ^ = -in,,  a ^ 2 ) = n3,,  — .2  , .(4), 


1 ^ = -sin  Xj  cos  X^  . i 
= (2/cosX)  sin  Xx , 


J>(1)  = - XDu> , <i>(2)  = D2u  = -X2W  + 

= -(2/cosX)  cos  Xx  , 


( i)  for  X^  ^ A . 


All  the  results  are  expressible  in  terms  of 

c^j  = (cos  Xi  cos  Xj)  1(cosXix  cos X jX > = (Xi  tanXi  - X^  tan  Aj)/(X2-X2). 


A double  integration  by  parts  of  cos  X.x  gives 


<w  . cos X . x>  = 
i 3 


,2 

1 *■*  A i 2 

— -K  (cosX.x  D‘u.)  = —S?  (a.COS  X.x>  + — (cosX.)c.  . 

X.  3 1 XT  1 3 X.  3 13 

3 3 3 


whence 


(cos A.)  <u>.  cos  X.x>  = - — s 2-  • 

3 1 3 x^-x^ 

i X j 

This  result  can  now  be  used  to  evaluate 

(♦j15  *j1)>  = XiXj  (Dwi  Do..)  = -XiXj(wi  D2(i)j  > = 

X.X3<u.(o.)  + 2X.A.(cosA.)  3 ( u.COSX  .x> 

1 3 1 3 13  3 1 3 

= X.X3  <w.w.>  - 4 X . X . c../(A2  -X2) 

13  13  1 3 13  x 3 


X i X j + 4X^X.  c^j/(A.  - Xj)  by  symmetry. 


Therefore  < w . o .) 

1 3 


(Xi-Xj)' 


<4’|1)4>f1)>  = - 4X.X.(X2  + X2)  c.  ./(X2  - X2)2  . 

1 3 131  3 13  1 3 

Similar  steps  using  (A2)  and  (A3)  give  the  results 


<«<2) 

.(2) 
♦j  ‘ 

= -8  X?  X2  c.  ./(X?  -X2)2 

13  13  1 3 

(A5) 

<*<1} 

A (3). 
*3  > 

= - 4X.X.  c.  ./(X2  - X2) 

(A6 ) 

4*2  c../u2 


(A7) 


<(* 


(2)  = . 


i 3 

Also  by  direct  integration 


->;> 


-l3’ 

.(3). 
♦j  > 

= 4 - 4(tanX^  + tan*j)/(*^  + * ^ ) 

(A8) 

"I4’ 

,(4) 

= 4c... 

i] 

( A9 ) 

(ill 


Direct  integration  gives 


<cos  x^y>  = 0 


(cos*i>  1<o, 


cos* 


1/3 ( X . cos* . ) ^ 

l l 


,y>  = 1/2*2 


> = 


1/2*| 


(A10) 


Using  these  results  we  find 


.(1), 

♦i  > " 

1 . 1 r Xi  * 2 _ ( 

2 3 '■cos*  J 

<*i3) 

. (3) 

♦i  > ' 

4/cos2xl  , < *<4>  4>  54 } : 

f— • 

" — ' -H 

-e- 

, (3). 

♦i  > ~ 

<*<2)  »<4)>  = i. 

(2)  „(2) 


(All) 


The  expreissions  quoted  for 


(*.(a'e) 

l 


,U,6) 


3 


3 (A12) 

(A13) 

in  Section  3 are  obtained 


immediately  from  these  results  since 

1.(4) 

2*i 


m3) 


,<3» 


1.(1) 

2*i 


1.(2) 

2*i 


1.  (4) 
2*i 
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Appendix  B.  Norms  of  eigenfunctions 
The  norms  are  found  from  the  expressions  of  Appendix  A by  writing 


Imtx^  tan  X^) 

Then  c.j  2 ( an  x ^ (He  X J ’ 


From  (A4) 

II 

*(1) 

II2  = 

u-n 

V>>  = -4  | X . | 2 

1 1 1 

f 2 (Re 

X?)  C.T 

1 11 

*i 

(2iImXi) 2 (2ReXi) 2 

whence 

II 

. (1) 
♦i 

H=I 

ix  i ! 

' (ReX2)  Im  ( X i tan 

Xi>  1 

1/2 

1 V 

( Im  X ^ ) 3 (Re  X . 

i3  J 

Similarly 

from 

(A5) 

II 

4<2) 

II  1 

ix  r 

~ Im  ( X^  tan  X ^ ) 

1 1/2 

*1 

II  2 

IV 

(Im  X . ) 3 ( Re  Xi)3 

and  from  (A8)  and  (A9) 


Also  from  (A3) 


' 2 Im ( X^  tan  X^) 
(Tm  Xi)  (Re  X^ 

2 Im  (X,  tan  X^) 
(Im  Xi)  (Re  Xi) 


1/2 

-1/2 


Asymptotically,  using  (3.5),  we  find 

II  ♦ II  , II  ♦1[2>  II  ^ \ kn/(j  in4kTi ) 3/2 

II  *'3)  II  , II  ♦J4’  II  ' 2/(*n4kV/2 


and 


(2  ku ) _1  ( |in4kn)“3/2 


for  large  k (in  fact  these  are  close  approximations  for  k > 2) . 


Appendix  C ( 1 ) Contour  integral  to  evaluate  £ 


2 2 
sin  Aj-sin  A^ 


Let 


.2  2 . 2 
T ( A A,  ) COS  A . 

(jjk)  ^ K ^ 


2 2 
2 sin  z - 2 sin  A, 


f(z)  = 


2 2 

(z  -a£)  C ( z ) 


where 


C(z)  = z + sin  z cos  z 
i 


and  consider  fUT  l dz  along  a contour  as  indicated  passing  to  the 


•i“ 


right  of  the  origin  but  to  the  left  of  all  the  zeros  A^  of  C(z)  with 


positive  real  parts. 


Since  f(z)  is  an  odd  function  of  z , the 
contributions  to  the  integral  from  the  two 
halves  of  the  imaginary  axis  vanish  and 


i 1“ 

/ f ( z ) dz  = j (residue  of  f at  z = 0) 


1 A. 


The  integral  also  equals  minus  the  sum  of  the  residues  at  the  poles  in 


Re  z > 0 , i.e.  at  the  zeros  A ^ of  C(z) 

2 

Since  C'(\.)  = 2 cos  A.  , the  residues  at  the  points  z = A . f A, 
3 3 3 k 


give  the  terms  of  the  required  series 


2/2  1 

The  residue  at  z = 0 is  sin  A /A  = = — , and  that  at  z = A 

V k cos^ A. 


is  -1/cos  A^.  Combining,  we  have 


-31- 


f 


l 

j 

j/k 

In  terms  of  = ( A ^ tan  A ^ - A^  tan  A ^ ) / ( A ^ - A2)  , the  result  is 

l -3k  I 

j(j^k)cos2A.  2cos2A, 

1 k 

(2)  The  second  result 

V Cik  / 2 2 

— 5 2 5 *1/4  A,  cos  A, 

(A2-A?)cos2a.  / k k 

3 K D 

is  obtained  by  the  same  method  applied  to  the  function  (z2  - A2)-1  f(z)  . 


. 2 .2 
sin  A . - sin  A, 

— 3 *_  = J- 

(A.  - A2)  cos2A  . 2 cos2a. 

J * J K 


1 


DAS/db 


-32- 


TABLE  2.  COMPUTED  COEFFICIENTS 


A (Collocation)  Case  1 


(c  ) real 


N = 5 

(o  = .1912 


61063369+00 

71418810-01 

34231232-01 

26835435-01 

21851789-01 


(c  ) imag. 
n 

59,-01) 

.61392605+00 
.13028745-01 
. 16279706-01 
. 42052902-01 
.59572337-02 


N = 10 

(a  = .745534,-02) 


61197602+00 

68466181-01 

25967889-01 

90125573-02 

21360312-02 

11323487-02 

19263364-02 

34329129-02 

14023023-01 

80038874-02 


.61710115+00 

.15751234-01 

-.12011278-01 

-.90003107-02 

-.48307601-02 

-.10557592-02 

.32729996-02 

.67216359-02 

-.15237260-01 

-.40843952-03 


N = 20 


(a  = .559897,-03) 


61143314+00 

69403359-01 

27184322-01 

10645565-01 

43879625-02 

17982206-02 

64708763-03 

11129095-03 

14400096-03 

26660636-03 

32713498-03 

36255633-03 

39570522-03 

44519822-03 

52774147-03 

64102925-03 

65523510-03 

55776629-04 

13495410-02 

38477966-03 


61566541+00 

14323567-01 

13490598-01 

10669075-01 

70334242-02 

45784616-02 

30236889-02 

20325576-02 

13870711-02 

95678757-03 

66326487-03 

45761730-03 

30660446-03 

18073945-03 

39125830-04 

19997518-03 

68723676-03 

13202449-02 

77284494-03 

17611918-03 


N = 40 

(a  = .470479,-03) 


.61140499+00 

-.69447217-01 

-.27231731-01 

-.10694934-01 

-.44391638-02 

-.18516945-02 

-.70369482-03 

-.17236297-03 

.76410530-04 

.18931400-03 

.23496628-03 

.24703267-03 

.24252237-03 

.23038100-03 

.21524926-03 

.19965022-03 

.18496796-03 

.17197961-03 

.16106763-03 

.15244672-03 

.14634454-03 

.14291931-03 

.14246104-03 

.14538541-03 

.15208382-03 

.16302043-03 

.17876976-03 

.19964429-03 

.22533701-03 

.25400598-03 

.28001656-03 

.28974349-03 

.25280907-03 

.10628609-03 

-.26422401-03 

-.99298726-03 

-.18573515-02 

-.78824002-03 

.24485574-02 

-.39475704-03 


.61558381+00 
.14243374-01 
-.13570348-01 
-.10748246-01 
-.71118366-02 
-.46559352-02 
-.31000388-02 
-.21075600-02 
-.14605170-02 
-.10285257-02 
- . 73 j60839-03 
-.52814852-03 
-.38252730-03 
-.27772044-03 
-.20135130-03 
-.14516082-03 
-.10350489-03 
-.72389746-04 
-.48971940-04 
-.31242143-04 
-.17615650-04 
-.68936447-05 
.19386252-05 
.99778445-05 
.18407256-04 
.28586092-04 
.42484144-04 
.62930762-04 
.94023941-04 
.14192214-03 
.21530492-03 
.32537540-03 
.48271942-03 
.68294941-03 
.85689188-03 
.72895106-03 
-.43119921-03 
-.29922575-02 
.64005520-03 
.39739060-03 


B (Fourier  Coefficients)  Case  1 


N = 5 

(a  = .124842 

.61196943+00 

.68648842-01 

.26453999-01 

.99103196-02 

.32725613-02 


N = 10 

(a  = .394030, 

.61154232+00 
.69252766-01 
.27044944-01 
.10524048-01 
.42858740-02 
.17166301-02 
.58801331-03 
.79233923-04 
.13855134-03 
. 195373j.3-0  3 


N = 20 

(o  = .619719, 

.61141044+00 

.69440008-01 

.27225579-01 

.10690154-01 

.44358049-02 

.18497402-02 

.70312456-03 

.17320371-03 

.74122276-04 

.18557761-03 

.22967359-03 

.23996520-03 

.23344906-03 

.21891356-03 

.20090414-03 

.18164388-03 

.16190454-03 

.14100743-03 

.11492843-03 

.65443454-04 


N = 40 

,-01)  (o  = .311031,-03) 


61740379+00 

.61139631+00 

.61555892+00 

16008436-01 

-.69460719-01 

.14218871-01 

11834770-01 

-.27246382-01 

-.13594793-01 

90707846-02 

-.10710136-01 

-.10772619-01 

55189619-02 

-.44547421-02 

-.71361037-02 

-.18675741-02 

-.46800459-02 

-.71985317-03 

-.31239312-02 

-.18884488-03 

-.21311882-02 

. 59547002-04 

-.14838687-02 

.17205230-03 

-.10515030-02 

.21719621-03 

-.75609051-03 

-02) 

.22855693-03 

-.55009736-03 

.22316833-03 

-.40385635-03 

61602638+00 

.20988085-03 

-.29839940-03 

14674471-01 

.19335072-03 

-.22131338-03 

13140750-01 

.17601934-03 

-.16430600-03 

10317183-01 

.15914360-03 

-.12172110-03 

66758288-02 

.14332582-03 

-.89639736-04 

42103435-02 

.12881870-03 

-.65302295-04 

26383389-02 

.11568607-03 

-.46732710-04 

16199345-02 

.10388763-03 

-.32501460-04 

93073679-03 

.93335177-04 

-.21562100-04 

42479420-03 

.83921507-04 

-.13138785-04 

.75533786-04 

-.66488299-05 

.68059953-04 

-.16534424-05 

.61396804-04 

.21797728-05 

-03) 

.55452139-04 

.50140227-04 

.51096406-05 

.73344634-05 

61560272+00 

14261664-01 

13552206-01 

10730131-01 

70936223-02 

46374686-02 

30811398-02 

20880438-02 

14402050-02 

10071175-02 

71073227-03 

50344894-03 

.45386793-04 

.41123137-04 

.37290717-04 

.33834981-04 

.30704110-04 

.27848840-04 

.25215059-04 

.22727375-04 

.20258110-04 

.17504634-04 

.13497920-04 

.46135928-05 

.90118913-05 
. 10264347-04 
.11192245-04 
.11879548-04 
.12397569-04 
.12812735-04 
. 13197701-04 
.13638738-04 
.14261768-04 
.15274646-04 
.16939768-04 
.17642350-04 

35550524-03 

24779171-03 

16767390-03 

10649439-03 

57950053-04 

16897873-04 

21634687-04 

57553125-04 


-34- 


_ . «.  • 

C (Biorthogonal  weighting  functions)  Case  1 


N = 5 

N 

= 40 

(a  = .101530, 

-01) 

(a  = • 

353545,00) 

61135616+00 

.61506104+00 

.61138740+00 

.61777331+00 

66043414-01 

.22716303-01 

-.70821297-01 

.72895306-02 

42069022-01 

-.30551851-01 

-.85014262-02 

.21413861-01 

44612018-02 

-.14961896-01 

-.53730856-01 

-.31473638-01 

52242245-02 

. 33942899-04 

-.41468516-01 

-.79426311-01 

.67799826-01 

-.11353737+00 

.26504809+00 

.44509070+00 

-.41282712+00 

-.22819515+00 

.79425870-01 

-.97633028-01 

N = 10 

.62999298-01 

.71562967-01 

(a  = .234948, 

00) 

-.65732111-01 

.14099306+00 

-.11796643+00 

.11325470+00 

61139172+00 

.61534385+00 

-.67857298-01 

-.16853370-02 

67526625-01 

.18427741-01 

.27470010-02 

.56969925-01 

36762522-01 

-.62455493-02 

-.78382837-01 

-.10058266+00 

53164362-01 

-.13168849+00 

.49289715-01 

.16432132+00 

18378647+00 

.13980365+00 

-.26295310+00 

-.31843493+00 

30478820-01 

. 30754357+00 

.27991951+00 

-.72133556-01 

49385241+00 

-.67001963+00 

.15114690+00 

.28860681+00 

55974969+00 

.32367657+00 

-.21247444+00 

-.14503381+00 

18598923+00 

-.20900748-02 

.15410253+00 

.59546073-01 

11616366-01 

-.12764101-01 

-.21342719-01 

-.23445248-01 

.85220357-01 

.11968018+00 

-.59676480-01 

-.34035228-01 

.95739569-01 

.16831081+00 

-.24286525+00 

.16911579+00 

-.22744965+00 

-.69862766+00 

N = 20 

.67420488+00 

-.73232653+00 

.76496596+00 

-.32548673+00 

(o=  .173240  , 

00) 

.31215408-01 

. 38959813+00 

-.80186294+00 

-.15218450+01 

61139365+00 

.61559773+00 

.16166298+01 

.22461629-01 

69515024-01 

.12498844-01 

.50750915+00 

.11626197+01 

23046318-01 

.45942632-02 

-. 51709488+00 

.22315959+00 

44470041-01 

-.82477171-01 

-.31367386+00 

.78136203+00 

93779008-01 

. 92207912-01 

-.12208629+01 

-.14669324+01 

76023699-01 

.83889777-01 

.11750934+01 

.27751912-01 

20072066+00 

-.42793861+00 

-.29149578-01 

.31095753+00 

46102554+00 

.41931194+00 

-.11234190+00 

-.42669927-02 

30125758+00 

-.83994967-01 

.56994109-02 

-.14144059-01 

29581786-01 

-.45013340-01 

18476482-01 

-.27760241-02 

73103741-02 

.50718615-01 

80956984-01 

-.47188486-01 

11736559-02 

-.10843889+00 

12660224+00 

.43565862-01 

18472532-01 

.11413889+00 

95335673-01 

-.39114366-01 

16368511-01 

-.19683665-01 

77247485-02 

.28738946-02 

73226797-03 

.93369492-03 

D (Optimal  Weighting  Functions)  Case  1 


N = 5 

N=  40 

(o  = .169764 

, -01) 

(o  = .535508, 

-03) 

61159911+00 

.61575004+00 

.61139783+00 

.61556352+00 

69100317-01 

.14443838-01 

-.69458680-01 

.14223256-01 

26839107-01 

-.13360062-01 

-.27244363-01 

-.13590491-01 

10296426-01 

-.10534172-01 

-.10708185-01 

-.10768378-01 

40524908-02 

- .68974781-02 

-.44528639-02 

-.71319176-02 

-.18658388-02 

-.46758996-02 

-.71817554-03 

-.31198313-02 

N = 10 

-.18726647-03 

-.21271293-02 

(o  = .775839 

,-02) 

.61069416-04 

.17349051-03 

- . 14798469-02 
-.10475154-02 

61141186+00 

.61557733+00 

.21852723-03 

-.75214178-03 

69429025-01 

.14239222-01 

.22984544-03 

- . 54618507-03 

27206685-01 

-.13575782-01 

.22440858-03 

- . 39998284-03 

10666254-01 

-.10755842-01 

.21104832-03 

- .29456071-03 

44089450-02 

-.71218203-02 

.19446362-03 

-.21750909-03 

18213547-02 

-.46682608-02 

.17710359-03 

-.16053709-03 

67402825-03 

-.31145654-02 

.16021549-03 

-.11799133-03 

14402424-03 

-.21241040-02 

. 14431612-03 

-.85945279-04 

10305702-03 

-.14789114-02 

.12977694-03 

-.61645959-04 

21401887-03 

-.10485142-02 

.11664314-03 

-.43116855-04 

.10474977-03 

-.28927634-04 

.94194293-04 

-.18034100-04 

.84737332-04 

-. 96600355-05 

.76310987-04 

-.32227700-05 

N = 20 

.68883314-04 

.17106726-05 

.62184615-04 

. 54778403-05 

(a  = .184497, 

-02) 

.56214307-04 

.83321986-05 

.50896950-04 

. 10469243-04 

61139748+00 

.61556317+00 

.46140892-04 

.12044575-04 

69459264-01 

.14222693-01 

.41847435-04 

.13176434-04 

27245032-01 

-.13591299-01 

.38072247-04 

.13956383-04 

10708896-01 

-.10769471-01 

. 34622139-04 

.14463590-04 

44535988-02 

-.71333188-02 

.31539193-04 

.14750772-04 

18665912-02 

-.46776211-02 

.28744154-04 

.14868894-04 

71894261-03 

-.31218762-02 

.26228844-04 

.14850977-04 

18804684-03 

-.21294947-02 

.23981761-04 

.14728459-04 

60276574-04 

-.14825254-02 

.21921008-04 

. 14523604-04 

17268593-03 

-.10504967-02 

.20033427-04 

.14259437-04 

21771170-03 

-.75541359-03 

.18366402-04 

. 13946097-04 

22901988-03 

-.54973397-03 

.16807824-04 

.13601162-04 

22357400-03 

-.40379472-03 

21020583-03 

-.,29862122-03 

19361435-03 

-.22180386-03 

17624872-03 

-.16505200-03 

15935618-03 

-.12271264-03 

14345350-03 

-.90859679-04 

12891212-03 

-.66740681-04 

11577716-03 

-.48379693-04 

-36- 


Table  2,  continued 
A (Collocation)  Case  2 


N = 5 N - 40 


(o 

= .377507,00) 

(o 

= .110522,00) 

. 607089184-00 

.66304201+00 

.58819868+00 

.63057330+00 

-. 35188187+00 

-.17427991+00 

-.36379328+00 

-.19661377+00 

- . 281575894-00 

-.47760766+00 

-.30736237+00 

-.46758012+00 

.473001774-00 

.40033579+00 

. 32520469+00 

. 37528259-01 

19663317400 

-.82746527-01 

.25972973+00 

. 38773020+00 

-.27555545+00 

. 94992820-02 

-.21074082+00 

-.33426678+00 

.25949360+00 

-.19895160-01 

. 18852364+00 

.30828986+00 

-.23995054+00 

.33259461-01 

-.16795803+00 

- .28014245+00 

N = 10 

.22568741+00 

-.35661868-01 

(o 

= .219051,00) 

.15300309+00 

-.21134837+00 

.26112026+00 

.40168331-01 

. 57906948+00 

. 61243860+00 

-.13950907+00 

-.24058998+00 

-. 36669589+00 

-.20909621+00 

. 19807889+00 

-.40232635-01 

30052660+00 

-.45171404+00 

.12780728+00 

.22368137+00 

. 27748704+00 

.20345909-01 

-.18469279+00 

.41250500-01 

. 19129017+00 

.29230234+00 

-.11673459+00 

-.20569544+00 

-.  23842'-74  + 00 

-.28631552-01 

.17154585+00 

-.39967122-01 

-.16859768+00 

-.26586103+00 

. 10660226  + 00 

. 18927678+00 

. 15557469+00 

-.11997603+00 

-.15770189+00 

. 39199884-01 

.23094832+00 

.23250230+00 

-.96262812-01 

-.17173319+00 

-.11012884+00 

.13362892-01 

.14413066+00 

-.36781278-01 

.86907761-01 

. 15505973+00 

-.12904927+00 

. 35049330-01 

-.76430163-01 

-.13659206+00 

N = 20 

.11462381+00 

-.31208296-01 

( o 

= .161249,00) 

.67303174-01 

-.97930284-01 

.11954374+00 

.29460705-01 

.58713482+00 

.62885981+00 

-.56595991-01 

-.98454383-01 

-.36300297+00 

-.19720814+00 

.81308153-01 

-.23017876-01 

-.30410254+00 

-.46249637+00 

.45556638-01 

. 81798684-01 

. 31862022+00 

.36416865-01 

-.65718552-01 

.22511513-01 

.25062907+00 

. 37256076+00 

-.39810450-01 

-.57370233-01 

-.26245858+00 

.83564809-02 

.37709177-01 

-.17224680-01 

-.19685610+00 

-.30897051+00 

. 16968374-01 

.27095119-01 

.23409013+00 

-.18709555-01 

-.19062765-01 

-.83537380-02 

.16493079+00 

.26660706+00 

.96351358-02 

.84214337-03 

-.20342026+00 

. 27805010-01 

-.18104902-02 

.14124047-02 

-.13808390+00 

-.22251441+00 

.16880088+00 

-.29508198-01 

. 10778646+00 

.17803678+00 

-.13921401+00 

.24179758-01 

Note:  If  the 

8 points  from  .465 

-.86363585-01 
. 98589820-01 

-.13958594+00 

-.30298234-01 

to  .535  are  omitted  in  calculating 

.63020841-01 

.77197284-01 

a , the  values 

obtained  are 

-. 39639047-01 
-.11717444-01 

.15390686-01 

-.85468226-02 

N = 5 o=  . 

298667 ,00 

.12554324-02 

-.23836171-02 

10 

990740,-01 

20 

210252,-01 

40 

792599,-02 

-37- 


B (Fourier  Coefficients)  Case  2 


N = 5 N = 40 


(0  = .736045,00) 

(0  = . 

218429, 00) 

. 56995639+00 

. 55534087+00 

. 58975801+00 

.63482024+00 

- . 38407921+00 

- . 2 7106706  + 00 

-. 36226089+00 

- . 192817  31+00 

31845103+00 

- . 55005451+00 

- . 30660976+00 

- . 46568250  + 00 

. 33589948+00 

-.  59793413-01 

. 32917903+00 

.41474951-01 

. 33028603+00 

. 25814066+00 

. 26445764  + 00 

. 396  30662+00 

-. 27822699+00 

. 1 34415  30-01 

- . 21375274+00 

-.  3389634  3+00 

N = 

10 

. 269  1 7150  + 00 

- . 16722648-01 

(0  = .207997,00) 

. 19744132+00 
- . 250501  36  + 00 

. 3252^495+00 
. 38513344-01 

. 58806410+00 

.62859458+00 

- . 17648978  + 00 

-.  29493297  + 00 

- . 36435445+00 

- . 19886986  + 00 

. 24432939+00 

-.34289289-01 

30825322+00 

- . 4 7186992+00 

. 167699  39  + 00 

. 28993511+00 

. 32808891+00 

. 34885158-01 

- . 2 3296141+00 

. 47865605-01 

. 26403508+00 

. 38890459+00 

- . 15  510928  + 00 

- . 26874551  + 00 

27766989+00 

. 45  313066-02 

. 22882351+00 

- . 41664174-01 

- . 21135522+00 

- . 35069076  + 00 

. 14980730  + 00 

. 267728  3 3 + 00 

. 27603655+00 

- . 34000815-01 

- . 22051420  + 00 

. 52543725-01 

. 21881788+00 

. 29652846+00 

- . 14094901+00 

- . 25047254+00 

- . 16217948+00 

-.10051733-02 

. 21755645+00 

-. 45204845-01 

. 13744334+00 

. 25206  361+00 

- . 21  10246  3 + 00 

.5531 3495-01 

N = 

20 

- . 13074123  + 00 

-. 23652538+00 

(0  = .191882,00) 

. 20871  39  3 + 00 

- . 46846117-01 

. 128175  30+00 

. 24036727  + 00 

. 58957268+00 

.63449720+00 

-. 20346676+00 

. 57303982-01 

- . 36264388+00 

- . 19315642+00 

- . 122995  36  + 00 

- . 22509122+0 0 

- . 3071 3315+00 

- . 46602409  + 00 

. 20139518  + 00 

-.  47197168-01 

. 32853490+00 

.41145275-01 

.12078234+00 

. 23157838  + 00 

. 26370712+00 

. 39600671  + 00 

-.  19735243  + 00 

. 59317960-01 

-. 27907255+00 

. 13194183-01 

-.  11704624  + 00 

- . 21480  368  + 00 

-. 21468521+00 

-.  33912911  + 00 

. 19490485  + 00 

- . 46003304-01 

. 26815637+00 

-.  16769374-01 

. 11428336  + 00 

.22585464+00 

.19634257+00 

. 32537525+00 

- . 19290000  + 00 

.62894641-01 

-.25169142+00 

. 38864  315-01 

-.  11341144  + 00 

-.  20314605  + 00 

- . 17778908+00 

- . 29426499  + 00 

. 18710000+00 

-. 40337969-01 

. 24288676+00 

-. 33183134-01 

. 10449237  + 00 

. 22795691+00 

. 16605129+00 

. 29165625  + 0 0 

- . 19881080  + 00 

. 77279381-01 

- . 23492995+00 

.50471317-01 

- . 1 3695378+00 

- . 174  314  33  + 00 

-. 15762143+00 

- . 26482061+00 

.92258291-01 

- . 264216  36-01 

.22529559+00 

- . 35674692-01 

. 14413953+00 

. 27715714+00 

- . 23148645+00 

. 68114654-01 

- . 16844438+00 

- . 22385806  + 00 

. 12317142+00 

-.  16466878-01 

-38- 


C (Biorthogonal  Weighting  Functions)  Case  2 


N - 

5 

N = 40 

(0  = . 

296387,00) 

(0  = . 

665455,01) 

58652741+00 

. 624  35478  + 00 

.58848195+00 

. 6 3296056  + 00 

33741776+00 

- . 16508999  + 00 

- . 38304045+00 

- . 24097565  + 00 

29684963+00 

- . 4 3551654  + 00 

- . 19047652+00 

-.  541  53776  + 00 

294  03796+00 

. 14827842+00 

. 75542744+00 

.16643047+01 

37653250-01 

-.10846512-01 

- . 24979798+01 

- . 285809  35+01 

. 22185785+01 

- . 14801681  + 01 

.39463186+01 

. 8912  3588  + 01 

N = 

10 

- .94  3562  34  + 01 

- . 10759307+02 

(0  = . 

307518,01) 

. 10056579  + 02 
- . 41267305  + 01 

.65261164+01 
. 14796249  + 01 

5886617  3 + 00 

. 63242231+00 

- . 26940930  + 01 

- . 18084892+01 

38162056+00 

- . 201001  33  + 00 

- . 56665776+00 

-. 37974272+01 

32037079+00 

- . 87871221+00 

. 40055308+01 

. 89648982+00 

13286031+01 

. 21618090  + 01 

. 12530144+01 

.45549932+01 

25776729+01 

- . 17168784+01 

- . 44096  398  + 01 

- . 95201063-01 

12810420+00 

-.4011411 3+01 

- . 24406790+01 

-.41197091+01 

60157050+01 

. 81086854+01 

. 18242068+01 

- . 10525019+01 

69125748+01 

- . 4 304665  3 + 01 

. 10884015  + 01 

- . 22752192  + 01 

26601005+01 

. 28426049  + 00 

. 37482455+01 

.41324295+01 

26820242+00 

. 17901652+00 

- . 24000814  + 01 

-.17623820+01 

. 24400731  + 01 

. 37477357  + 01 

- . 36617273  + 01 

- . 36  56104  3 + 0 1 

N = 

20 

. 40659974  + 01 

. 17940227  + 01 

(a  = . 

240200,01) 

- . 56689912  + 00 
.82535621+00 

. 79464592  + 00 
. 21  386242  + 01 

58854527+00 

.63018823+00 

- . 29357848  + 01 

. 3 3170  36  3 + 01 

362542  39  + 00 

- . 16  309715  + 00 

- . 421  18856  + 01 

- . 97  371110  + 01 

39442898+00 

- . 79456420  + 00 

. 89781090  + 01 

. 10282159+02 

90945917+00 

. 11895634  + 01 

- . 10  300884  + 02 

- . 50393864  + 01 

12321781+01 

- . 10696868+01 

. 33223339+00 

. 59910794  + 01 

75953955 +00 

-.  11929522  + 01 

-.  12261743  + 02 

-. 23590204+02 

24599904+01 

. 53198091+01 

. 25161335  + 02 

. 33406487  + 00 

58505628+01 

- . 57619814  + 01 

. 79291270  + 01 

. 18463717  + 02 

44671741+01 

.18077749+01 

- . 84301004+01 

.29131 306+01 

95020428+00 

. 51546749  + 00 

- . 44095077  + 01 

. 12657680  + 02 

46288896+00 

- .66689474-01 

-.  19545304  + 02 

-.  23355228  + 02 

33980731+00 

-. 14919883+01 

.18855635+02 

.60294698+00 

19370439+01 

.97540758+00 

- . 59048611  + 00 

. 49891589+01 

31444669-01 

. 19517805  + 01 

-.18217431+01 

- . 10922078+00 

21437625+01 

- . 721  33705  + 00 

. 10474451  + 00 

- . 2 32  36669  + 00 

26357984+00 

- . 19435165  + 01 

17111665+01 

. 81203042  + 00 

37875147+00 

. 39581158  + 00 

21841977+00 

- . 12684664  + 00 

49799360-01 

-.27011399-01 

-39- 


D (Optimal  Weighting  Functions)  Case  2 


N 

= 5 

N 

= 40 

(0  = , 

.121  3 8 2, 01  ) 

(O  = .204838,01) 

59108404+00 

. 63404593+00 

. 588  3 3995  + 00 

.63071703+00 

35984104+00 

- . 19  119  311  + 00 

- . 36445971+00 

- . 19688687+00 

30363506+00 

- . 46219734  + 00 

- . 30895874+00 

- . 469762  36  + 00 

33234702+00 

. 46  340538-01 

. 32678397+00 

. 37388938-01 

26760696+00 

. 40214169+00 

. 26205  347  + 00 

. 39221676+00 

- . 28062099  + 00 

.9 3487867-02 

- . 21612499  + 00 

-.3430591 7+00 

N 

= 10 

. 26682878+00 

- . 20821836-01 

(0  = . 

862277, 00) 

.1951 3340+00 
- . 25276970  + 00 

.32115053+00 
. 34401735-01 

58983302+00 

. 6 32  17881  + 00 

- . 17871569  + 00 

- . 29905374+00 

36149827+00 

- . 19510081+00 

. 2421481  1 + 00 

-.38422405-01 

30532441+00 

- . 46791507  + 00 

. 16556391+00 

. 28578571+00 

33073473+00 

.39232946-01 

- . 23504967+00 

. 4 369492  3-01 

26612317+00 

. 39403208  + 00 

- . 15714877+00 

- . 27294209+00 

27654766+00 

. 11  122597-01 

. 2268335  3 + 00 

- . 45892067-01 

21211596+00 

- . 34  1 3 3417  + 00 

. 14786680  + 00 

.26346264+00 

27073417+00 

-.  19149920-01 

- . 22240446  + 00 

.48231564-01 

19891 346+00 

. 32276706+00 

- . 14278787  + 00 

- . 25483952+00 

24912561+00 

. 35961911-01 

. 21576945+00 

-.49636553-01 

. 13570904+00 

. 24755440+00 

-.  21270471  + 00 

. 50711810-01 

N 

= 20 

-.1  32  36  3 38  + 00 

- . 24123280+00 

(0  = . 

, 1 33466,01  ) 

. 20714962+00 
. 12667  339  + 00 

- . 51680766-01 
.23538134+00 

58805975+00 

. 63042476+00 

-.20489849+00 

. 52138400-01 

36506118+00 

- . 19726869  + 00 

-.12435064+00 

- . 2 3046878  + 00 

30974839+00 

- . 47016876  + 00 

. 20012298  + 00 

-. 52834260-01 

32587561+00 

. 36979368-01 

. 11960960+00 

. 22562209+00 

26107277+00 

. 39181367  + 01 

- . 19840117  + 00 

. 52968735-01 

28164125+00 

. 89569794-02 

- . 1179  388  3 + 00 

- . 22164478+00 

21716120+00 

- . 34343723  + 00 

. 19421522  + 00 

-.53480232-01 

26579367+00 

-.  21184948-01 

. 11388429+00 

.21753837+00 

19411167+00 

. 32080286+00 

- . 19285919  + 00 

.53435815-01 

25376931+00 

. 34069547-01 

- . 11265369  + 00 

- . 21421830+00 

17968707+00 

- . 299  37065  + 00 

. 1891  3580  + 00 

- . 5 3816912-01 

24120909+00 

- . 38724  386-01 

. 10911814+00 

. 21067804+00 

16465991+00 

. 28549818  + 00 

- . 18804207+00 

. 53672894-01 

23591712+00 

.43421326-01 

- . 10819491  + 00 

- . 20784126+00 

15797897+00 

- . 27  320231+00 

. 18469257  + 00 

- . 5 395  782  3-01 

22604069+00 

-.46139496-01 

14711094+00 

. 26322744  + 00 

22312404+00 

. 48008004-01 

NOTE:  Omitting 

the  8 points 

14347210+00 

- . 25505199+00 

between  .465 

and  .53  gives 

21511944+00 

- . 49838479-01 

N = 5 

0 = .116627,01 

10 

.842872,00 

20 

.132231,01 

40 

. 203624 ,01 

-40- 


RHS  for  Case  D 2 


( cl  ) , < (d  ). 

n real  n imag 


. 527126  34+00 

.44575244+00 

- . 42918404  + 00 

- . 34640071+00 

- . 3 ‘■'688970  + 00 

- . 58352952  + 00 

. 30454501+00 

-.  51253072-01 

. 24841781+00 

. 32093542+00 

- . 2893655  3 + 00 

-.49650438-01 

22194791+00 

- . 393091  37  + 00 

. 26283596+00 

- . 64103497-01 

. 19233676  + 00 

. 28  309255  + 00 

- . 25475471+00 

. 48408798-03 

19013182+00 

- . 32962089+00 

. 24114218+00 

- . 66229125-01 

. 16486112+00 

.26028875+00 

- . 23552364+00 

. 20156783-01 

-.15744683+00 

- . 294  799  32  + 00 

.22667289+00 

-. 66291870-01 

. 14781493+00 

. 24433808  + 00 

-. 22236924+00 

. 302  31793-01 

- . 14268221+00 

-. 27183997+00 

. 21593248+00 

-.65743316-01 

. 13591914+00 

. 23225148  + 00 

- . 21245583+00 

. 361  35656-01 

- . 13208243+00 

- . 25514869  + 00 

. 20745715+00 

- . 64994085-01 

.12700295+00 

. 22262010  + 00 

- . 20455069+00 

. 39884693-01 

- . 12398774+00 

- . 24225424  + 00 

. 20049833+00 

- .64186285-01 

. 11999510+00 

. 21467241  + 00 

- . 19800745+00 

. 42  393384-01 

- . 11753854+00 

-.  23187091  + 00 

. 19462068+00 

- .63379679-01 

. 11429367+00 

. 20794512  + 00 

- . 19244695+00 

. 441  30015-01 

- . 11223951+00 

-.22325366+00 

.18955116+00 

-.62597247-01 

. 10953394+00 

. 2021  3.854  + 00 

- . 18762645+00 

.45361058-01 

- . 10777996+00 

- . 21593749  + 00 

. 18510641+00 

-.61849424-01 

TABLE  3 


Computed 

Coefficients  for 

Zero-End  Displacement 

Problem  (E4) 

(c  ) , 

n real 

(c  ) . 
n imag 

N = 

5 

N 

= 40 

. 11490602-01 

. 16880880  + 00 

.84393964-02 

. 16735048  + 00 

- . 16110192-01 

.66001308-01 

- . 19706981-01 

. 65  384974-01 

- . 14339124-01 

. 41379458-01 

- . 18124259-01 

. 41190661-01 

-.11674771-01 

. 30448630-01 

- . 15481876-01 

.30540044-01 

- .95871504-02 

.24290585-01 

-.13331680-01 

. 24582611-01 

- . 11667422-01 

. 20769061-01 

- . 10364452-01 

. 18107270-01 

N = 

10 

- . 93226098-02 

. 161  34865-01 

. 96778788-02 
- . 1817083  3-01 
-.  16417219-01 
-.13676391-01 
- . 11474775-01 
- .97904070-02 
-.84885197-02 
- . 74622143-02 
-.66370235-02 

- . 59616720-02 

. 16817558+00 
. 65860470-01 

. 41481841-01 
. 30699111-01 
. 246  38  308-0  1 
. 20740474-01 
. 18008303-01 
. 15976204-01 
. 14398375-01 
. 13132877-01 

- . 84720298-02 
-. 77646417-02 
- . 71668816-02 
-.66548045-02 
- .621094  32-02 
- . 58222854-02 
-. 54789477-02 
-.51732931-02 
- . 48993267-02 
-.46522762-02 
- . 44282941-02 

. 14608224-01 
. 1 3 386980-01 
. 12384557-01 
. 1 1544617-01 
. 10828878-01 
. 10210389-01 
. 96696002-02 

. 91919710-02 

. 87664474-02 
. 8 3844685-0  2 
. 80392931-02 

- . 42242437-02 

. 77255377-02 

N = 

20 

- . 40375399-02 

. 74388442-02 

-. 38660330-02 

. 71  756468-02 

.88332921-02 

. 16768864+00 

- . 37079195-02 

. 69329957-02 

-. 19208893-01 

.65608801-01 

- . 35616747-02 

.67084283-02 

- . 17557180-01 

. 41 356522-01 

- . 34260004-02 

.64998724-02 

-. 14865636-01 

. 30663542-01 

- . 32997837-02 

.63055693-02 

- . 12680160-01 

. 24671699-01 

-. 31820646-02 

.61240146-02 

-. 10990873-01 

. 20828860-01 

- . 30720113-02 

. 59539179-02 

- . 96706884-02 

.18141507-01 

- . 29688973-02 

. 5794  1600-02 

-. 86177045-02 

.16146455-01 

-.28720870-02 

. 56437680-02 

- . 77607587-02 

. 14599547-01 

-.27810212-02 

. 55018915-02 

- . 70507870-02 

. 1 3 3600  38-01 

- . 26952048-02 

. 5 3677824-02 

-. 64534499-02 

. 12341069-01 

- . 26142000-02 

. 52407815-02 

- . 59441896-02 

. 11486079-01 

-.25376154-02 

. 51203039-02 

- . 55050490-02 

. 10756609-01 

- . 24651012-02 

. 50058289-02 

- . 51226218-02 

. 10125559-01 

- . 23963454-02 

.48968929-02 

- • 47867055-02 

.95732568-02 

- . 23310661-02 

. 47930786-02 

-.44894039-02 

. 90850566-02 

- . 22690091-02 

. 4694012  3-02 

- . 42245104-02 

. 86498135-02 

- . 22099449-02 

. 45993563-02 

- . 39870794-02 

. 82588878-02 

- . 21536654-02 

. 45088055-02 

- . 37731218-02 

. 79  .-4692-02 

- . 35793839-02 

. 75841130-02 

Papkovitch-Fadle  Functions 


Figure  2 (a) 


Papkovitch- Facile  Functions 
Im(tanX  cos  Ax  - x sin  X 


Figure  2(b) 


Figure  3 
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Figure  4 
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such  equations  that  have  been  proposed  in  the  past  are  unstable  with  respect 
to  the  order  of  truncations— — 

A systematic  search  among  all  possible  weighting  functions  in  the  bi- 
orthoqonal  family  for  use  in  a Galerkin  method  leads  to  a unique  choice  of 
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Johnson  and  Little  (1965),  and  with  a method  of  direct  collocation.  The 
method  of  optimal  weighting  functions  produces  considerably  more  stable  values 
for  coefficients  as  the  order  of  truncation  is  changed.  This  opens  the  way 
for  an  examination  of  the  convergence  of  the  eigenfunction  expansions  for 
data  of  this  class,  and  should  also  prove  useful  for  problems  involving 
cracks. 


